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Abstract 

A nearest neighbour spin pair of the quasi-two-dimensional three-state Potts 
model interacts with the strength J(> 0) in the a;?/-plane and with \J (0 < 
A -C 1) in the z-axis. The phase transition is of second-order when A = 
and is of first-order when A > 0. The dimensional crossover occurs with a 
change of the order of the phase transition. We study the finite-size effects 
of the phenomenon by using a Monte Carlo method with a multi-spin coding 
technique. The prediction of the finite-size scaling theory is consistent with 
the Monte Carlo results. 



1 Introduction 



A quasi-two-dimensional system is a three-dimensional one in which the ra- 
tio A of the interplanar to the intraplanar exchange interactions is small. If 
the system exhibits a phase transition, we can see a crossover from two-di- 
mensional to three-dimensional behaviour as the critical point is approached. 
Many authors have studied quasi-two-dimensional antiferromagnets by neu- 
tron scattering experiments f^. 

The theoretical studies of the dimensional crossover have been carried 
out by using perturbation theory high temperature series expansion 

[^]-[]T^, generalized homogeneous functions |T^, and rigorous approach ]TB|- 
[p^ . The results are that the critical and the crossover temperature are 
singular with respect to A. They behave as A^/"^ where is the crossover 
exponent It has been shown that is equal to the critical exponent 7 
for the susceptibility of the system with A = 0. The present author discussed 
finite-size scaling and performed Monte Carlo simulaltions of the quasi-two- 



dimensional Ising model for the first time [|19 



In the above mentioned researches there was an assumption that the order 
of the phase transition was unaltered. In this paper we study finite-size 
scaling of the dimensional crossover in which it changes. We perform Monte 
Carlo simulations of the three-state ferromagnetic Potts model pO|. The 



phase transition is of second-order in two dimensions and is of first-order 



in three dimensions |22]. In the next section we review the finite-size scaling 



theory for the quasi-two-dimensional systems briefly. The finite-size scaling 
form of an effective transition temperature is presented. In section |^ we 
describe an algorithm for the multi-spin coding technique used in our Monte 
Carlo simulations of the quasi-two-dimensional three-state Potts model. The 
Monte Carlo data are compared with the prediction of the finite-size scaling 
theory in section ^ A summary is given in section ^. 



2 Finite-size scaling theory 



We review the finite-size scaling theory for the quasi-two-dimensional sys- 
tems briefly. The detailed discussion is in the reference |19|. Let us consider 



the three-state Potts model on the simple cubic lattice for the sake of con- 
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creteness. The Hamiltonian is 



= E [1 - ^(«^- ^.O] + ^ E[i - ^(^- 1)] (1) 

where cxj is a Potts spin variable located iih lattice site and which takes on 
the value 1, 2, and 3. The first summation is over all nearest neighbour pairs 
on the lattice, the second summation over all lattice sites. The strength Jij 
of the interaction for the nearest neighbour pair ij is J(> 0) in the a;?/-plane 
and AJ (0 < A -C 1) in the z-axis. H is an external magnetic field. When 
A = 0, the equation (||) consists of two-dimensional three-state ferromagnetic 
Potts models which are independent of each other. 

We assume that the free energy per spin measured by k-oT , where k-Q is 
Boltzmann's constant and T is the temperature, is a generalized homoge- 
neous function of variables 

to = T/T(0) - 1, 

h = H/k^T, and A as Iq, h, X ^ 0, where T(0) is the critical temperature of 
the system with A = 0. We can derive 

f{to,h,X) = |top-"/^(V|io|^+^A/|to|^) (2) 

with a scaling function f^{x,y) where -|- (— ) refers to to > (to < 0), 
and a + 2/5 + 7 = 2 where a, /3, and 7 are the critical exponents of the 
two-dimensional system for the specific heat, the magnetization, and the 
susceptibility, respectively. The number is the crossover exponent and is 
equal to 7. From we get the behaviour of T(A) that is the transition 
temperature of the system with A as follows. 

r(A)/T(0)-l = A^A^/-^ (3) 

where is a constant. 

Let us consider the three-state Potts model (P on an L x L x L simple 



cubic lattice to see the finite-size effects |^-|]25| of (y). To avoid surface 
effects we impose periodic boundary conditions. We assume that the free 
energy per spin is a generalized homogeneous function of to? A, and L as 
to, A, 1/L ^ and the system is characterized by L/^ito) where ^(to) is 
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the correlation length of the system oih = X = l/L = 0. Using (Q) we can 
derive 

/(to, h, X, L) = /(to h a L*/^) (4) 

where f{x,y,z) is a scaling function and u is the critical exponent for ^(to). 

Let us define an effective transition temperature as the position, Tl{\), 
of the peak of the specific heat. From (^ we get 

ri(A)/T(0) - 1 = L-'/" T(A L-^/") (5) 

where T(x) is a scaling function. If T(x) x^^'^ as x — +oo, the 

equation (^ is reproduced in the hmit L +oo for a fixed value of A(> 0). 
The prediction (^) will be compared with the Monte Carlo data in section ^ 



3 Monte Carlo simulations 

To confirm the prediction (Kl) of the finite-size scaling theory, we perform 



Monte Carlo simulations p^, ^ of the quasi-two-dimensional three-state 
Potts model (|I|) in if = on the L x L x L simple cubic lattice under fully 
periodic boundary conditions. We use a multi-spin coding technique ||28|, ^ 
to simulate a large number of systems simultaneously. Since the FORTRAN 
compiler on the HITAC S-820/80 computer, we have used, treats 32-bit in- 
tegers, we can update 32 systems independently. Three-state Potts spin 
variables located at identical lattice sites are stored in the 32-bit positions of 



two words 30, 31 



An algorithm is as follows. Let us consider a flip of a spin ctq. The change 
in the energy on flipping the spin, ctq — ctq, is 

AE/J = -j2{[l-Siao,a,)]-[l-6{a'o,a,)]} 
-xj2{[l-S{ao,a,)]-[l-6ia'„a,)]} 

k=5 

where aj (j = 1,2,3,4) is the nearest neighbour spin of ao in the xy-plane 
and (Tfc {k = 5, 6) in the 2;-axis. Using variables 



n. 



xy 



^j:{[^-S{ao,a,)]-[l-S{aia,)] + l} 
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and 



we have 
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"-2 = ^{[1 ~ '^(^0, CTfe)] - [1 - S{(7q, (Tfc)] + 1}, 

-AE/J = n^y + Xuz- A- 2X. 



Defining a variable n = bn^y + riz + 10, we may see the energy change as a 
function of n: {-AE/J)n, n G {10, 11, . . . , 54}. 

Using a random variable W which takes an integer value, we flip the spin 
if n + > 32 where w is a possible value of W. The distribution of W is 
determined by 

Prob{iy > 32 - n} = min(e^(-^^/-^)", 1) 

r ^K{-AE/j)„^ 10<n<32, 
~ j 1, 32 < n < 54, 

where K = J/k-^T. Thus the update is accepted with the probability 
min(e~^'^^, 1) where P = 1/k^T. This is the same procedure as of the 



Metropolis algorithm [Q. We express n + w {10, 11, . . . , 76}), which 
can be calculated with logical operations, as the binary notation: I]f=o^'2', 
Xi G {0, 1}. In our algorithm we carry out the update when X5 = 1 or 
xq = 1 although the Metropolis procedure needs to refer to the inequal- 
ity r < min(e~^^^, 1) where r is a possible value of a random variable R 
with uniform distribution over [0,1]. The algorithm explained here is for 
A G (0,1/4). It is sufficient to study the system since we are interested in 
the case of the small value of A. 

The pseudorandom numbers are generated by the Tausworthe method 
[ ^ , |3^ . We measure physical quantities at a temperature over 10^ Monte 
Carlo steps per spin (MCS/spin) after discarding 10"^ MCS/spin to attain 
equilibrium. Physical quantities are calculated by the multi-step bitwise 



summation algorithm |3^, |3J]. Our algorithm achieves a speed of 25 million 
spins per second on a 30 x 30 x 30 lattice with measurements at every step. 
Let us denote the average of a physical quantity, O, in each system by 
i = 1, 2, . . . , 32. The expectation value is given by 

1 32 
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the standard deviation by 



A{0) = ({0Y-{0f)"" /V^. 



4 Monte Carlo results 



In the reference we reported that there was hysteresis in magnetic quanti- 
ties (magnetization, susceptibihty, . . . ) but we could not see it in the energy, 
the specific heat: 

and the fourth-order cumulant of the energy for the small value of A or a 
small number of MCS/spin. The hysteresis vanished for the large value of 
A or a large number of MCS/spin. We have observed similar behaviour for 
the quasi-two-dimensional three-state Potts model. According to [|T^], we 
analyse the peak position, Tl{X), of the specific heat since we can obtain 
useful information about phase transitions from it. 

Figure 1 shows the temperature dependence of C of L = 20 system for 
various A. The solid curves are obtained by the smoothing procedure of the 



fourth-order i?-spline |]35|. As A increases, the shape of the curve becomes 
sharper and the peak position T2o(A) becomes higher. Then we get the peak 
height and T2o(A). 

Figure 2 shows the finite-size scaling plot of Tl{X). We have used the 
exact values of the critical temperature, k-oTifi)/ J = l/ln(l -|- v^), and the 
critical exponents, = 7 = 13/9 and = 5/6, of the two-dimensional three- 
state ferromagnetic Potts model |2^. The data fall on a common curve. It 
is consistent with (|^). The slope seems to approach 1/0=1/7 = 9/13. It 
indicates that the asymptotic behaviour of the finite-size scaling function is 
T(x) const. X x^^'^ as X +00. 



5 Summary 

The dimensional crossover occurs with a change of the order of the phase 
transition in the quasi-two-dimensional three-state Potts model. We have 
studied the finite-size effects of the phenomenon by using a Monte Carlo 
method. An algorithm of a multi-spin coding technique for this model has 
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been presented. The prediction of the finite-size scahng theory has been 
consistent with the Monte Carlo results. We have seen the asymptotic be- 
haviour of the finite- size scaling function on the system with A < 0.24 and 
L < 30. 
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Figure captions 



Figure 1 Temperature {K — J/ k^T) dependence of the specific heat of the 
L = 20 system for various A: 0.04, 0.08, 0.1, 0.12, 0.14, 0.16, 0.18, 0.20, 
0.22, 0.24. The sohd curves are obtained by the smoothing procedure 
of the fourth-order i?-sphne. As A increases, the shape of the curve is 
sharper and T2o(A) becomes higher. 

Figure 2 Finite-size scahng plot of the effective transition temperature. The 
data of the system with L = 10, 20, and 30 are denoted by Q) ^-^d □, 
respectively. We have set that A;Br(0)/ J = l/ln(l+\/3), = 7 = 13/9, 
and u = 5/6. The data fall on a common curve. The solid line shows 
that (9/13) ln(AL^^/^^). Errors are less than the symbol size. 
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